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Abstract 

We compute the distributions of charged particles at large transverse momenta in 
pp(p), pA and AA collisions in the framework of perturbative QCD, by using collinear 
factorization and the modern PDFs and fragmentation functions. At the highest cms- 
energies the shape of the spectra measured in pp(p) collisions at large qx can be well 
explained. The difference between the data and the lowest-order computation is quan- 
tified in terms of a constant fCfactor for each energy. The i^-factor is found to sys- 
tematically decrease with growing ^Js. Also a lower limit for the partonic transverse 
momentum, po, is extracted for each yfs based on the comparison with the measure- 
ments. A systematic increase of po as a function of is found. Nuclear effects in the 
charged-particle spectra in pA and AA collisions at RHIC and LHC are studied in the 
framework of collinear factorization by applying the EKS98 nuclear corrections to the 
parton distributions. The nuclear effects are shown to mostly enhance the computed 
spectra. A comparison with the recent PHENIX data from central and peripheral 
Au+Au collisions at RHIC is done. 
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1 Introduction 



Perturbative QCD (pQCD) and collinear factorization has been shown to be successful 
in explaining the production of observable jets in pp collisions Jl], |2|. The elements 
in such calculations are the parton distribution functions (PDF) which contain non- 
perturbative input at some initial scale perturbatively computed cross sections for 
parton-parton scatterings, and the infrared-safe measurement functions which deter- 
mine the cross section to be computed. For the same energies it is, however, surprisingly 
difficult to predict inclusive one-particle distributions of charged hadrons in hadron- 
hadron collisions from perturbative QCD, even at large cms-energies i/i > 100 GeV 
and px ^> Aq CD , where the perturbation theory should work very well. 

In the framework of collinear factorization, the inclusive charged-particle transverse 
momentum spectra can be computed by convoluting the inclusive production of a 
parton k further with the corresponding fragmentation function of k into a hadron h, 
Dk^h- Schematically 

d(T PP-*h+x = J2fi(xi,Q 2 )fj(x2,Q 2 ) ®a i ^ k+x {x u x 2 ,Q 2 ,a s {ii 2 )) ®D^ h {z,^ F ), (1) 

ijk 

where the partonic cross sections a l i^ k+x are expressed as a power expansion in a s (fi). 
The factorization scale is denoted by Q, the renormalization scale by /x, and the frag- 
mentation scale by hf- State-of-the-art calculations where next-to- leading order (NLO) 
partonic cross sections and NLO PDF are incorporated with NLO fragmentation func- 
tions H |J extracted from e + e~ collisions can be found in [^, [J. While in these calcu- 
lations the NLO results are well compatible with the measured cross sections J7], [8|, ||] 
at large transverse momenta, the shape of the data is generally not well reproduced 
& t It ^ 5 GeV. This can be thought to be as a consequence of the fact that the frag- 
mentation functions of gluons are not sufficiently constrained by the e + e~ data, or 
that the effects of intrinsic &t both in the fragmentation functions and in the PDF 
are neglected. Especially at lower cms-energies the fey-effects are expected to become 
important, and phenomenological approaches studying this problem starting from pp 
collisions at y/s ~ 20 GeV have been developed e.g. in QTO], [T| based on the lowest-order 
(LO) pQCD approach. 

In this paper, we will not make an effort to perform a NLO computation of the 
charged-particle transverse momentum spectra, or to invoke a phenomenological model 
for the intrinsic fey-effects. Instead, we compute the inclusive production of charged 
particles in LO but strictly remaining within the collinearly factorized approach, sup- 
plementing the computation only with a constant i^-factor but no other phenomeno- 
logical factors. We show that the shape of the computed transverse momentum distri- 
butions agrees reasonable (surprisingly?) well with the measured ones in most cases. 
Emphasizing the highest values of transverse momenta measured at each cms-energy, 
we determine the ^-factors (K = da^ xp / da^ Q ) based on a \ 2 fit to the data. A sys- 
tematic decrease of K with growing cms-energy is found. 
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Particle production in hadronic and nuclear collisions is often modelled in terms of 
minijet production above some transverse momentum cut-off po, supplemented with a 
nonperturbative component for p? < Po [0> 0- The l° wer limit of partonic trans- 
verse momentum needs to be phenomenologically determined from the measurements. 
This is done e.g. using the total and inelastic cross sections in pp and pp scatterings 
in context with eikonal models |L2|, [15|]. The charged-particle spectra measured at 
midrapidities, however, offer another independent way of obtaining constraints for the 
cut-off scale. Our second goal in this paper is to try to find lower limits for p based 
on the comparison of the computed and measured px distributions of charged hadrons. 
A systematic increase of po as a function of ^/s is found. 

Our third goal is to study the magnitude of the nuclear effects in high-p^ hadron 
production in pA and AA collisions at RHIC and LHC energies by strictly keeping 
the framework of collinear factorization and DGLAP evolution. For doing this we 
apply the EKS98 nuclear modifications [|16] to the PDF, and investigate the relative 
difference between AA collisions and the corresponding pp collisions. Both RHIC and 
LHC energies are studied. We show that the pQCD spectra in central Au+Au collisions 
at RHIC at qr < 10 GeV are enhanced by the nuclear effects at most by 15%. 

The fourth and final goal of the present paper is to compare the computed pQCD 
spectra of charged hadrons to the very exciting recent data from RHIC measured 
by PHENIX [fT? |. To obtain the absolute normalization for the spectra, we need to 
extract the if-factor for y/s = 130 GeV based on the obtained behaviour K(y/s). Two 
centrality classes are studied. The more peripheral one is found to be consistent with 
pp collisions, while in the 0...10% central sample, the computed transverse momentum 
spectra systematically lie above the data - especially at the highest transverse momenta, 
where the emphasis of our approach is. 

The paper is organized as follows. In Sec. 2 we present the formulae for the computed 
pQCD cross sections in detail. Sec. 3 explains how the i^-factors are extracted and 
summarizes the pp(p) results. In Sec. 4 the study is extended to nuclear collisions. 
Sec. 5 contains the comparison with the PHENIX data. Discussion is given in Sec. 6. 



2 The formulae 

The formulae for computing inclusive charged-particle production in LO pQCD are 
given below. In LO, the partons are produced back-to-back in the transverse plane 
according to the differential cross section 

dcr AB^ki+x da ij ^ k K 

, 2 , , = J2 X lh/A( x l,Q 2 ) x 2fj/B(x2,Q 2 )-^ (S,t,u) (2) 

dp A T dy x dy 2 ij ' dt 

where A and B denote the colliding hadrons or nuclei. The rapidities of the final state 
partons k and I are labelled by y\ and y 2 , and transverse momentum of each parton by 
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Pt- The fractional momenta of the colliding partons i and j are x 1>2 = ffle** + e***), 
i.e. the incoming partons are collinear with the beams. The parton distributions are 
obtained from the CERN-PDFLIB library |TB[, and the factorization scale is Q ~ px- 
The Mandelstam variables for the subprocesses are denoted by s, t and u. In LO, 
the partonic cross sections are da/dt ~ a^(fi 2 ). The 1-loop expression for the strong 
coupling constant is used here, so 



a s (i/ 



12tt 



(33 - 2N f ) log^/A^, 



(3) 



where /i ~ px is the renormalization scale. In the present study all partons are consid- 
ered massless and the number of active flavours, Nf = 5,4,3 below the corresponding 
threshold scales Qq = m t = 174 GeV, Q$ = nib = 4.5 GeV and Q4 — m c = 1.5 GeV. 
For Nf = 4, we will use the Aq^ d as given by PDFLIB [18| for the chosen set of parton 



distributions. The change of Aq^^ is computed by requiring matching of a s at the 
threshold scales. The cross sections for the eight different partonic subprocesses can 
be found e.g. in [l]| . 

Inclusive cross section for production of a parton of a flavour / and a rapidity yf is 
obtained by integrating over one of the rapidities in Eq. (Q) and keeping track of the 
parton flavours [20], 



da AB->f+X 



dp T dyf 



da 



AB^kl+X 



(ki) 



X 



(ij)(kl) ■*• 1 ° kl •*• 1 "V 
X Q 2 ) X 2 fj/B{x2, Q 2 ) 

+%lfj/A(Xl, Q 2 ) X 2 fi/B(x 2 , Q 2 



kl 



da ij ^ k l da ij ^ kl ~ 

—z (s,t,u)5 fk + —j (s,u,t)S f i 
dt dt 

da ij r: k i . r da ij r H ~ ^' 
{s,u,t)df k + —j (s,t,u)dfi 
dt dt 



'(5) 



where the summations run over the pairs, i.e. (kl), (ij) = gg, gq, gq, qq, qq, qq, qq, (q = 
u, d,s, . . .) without a mutual change. 

For inclusive hadron production through fragmentation of the parton /, let us define 
the transverse momentum of a hadron h be qx, its transverse mass mx and rapidity y. 
We define the fraction z to be the ratio of the energy of the hadron and the energy of the 
parton / from which the hadron originates: z = Eh/Ef. The hadron is considered to 
be produced collinearly with the parton. With these definitions the partonic variables 
become related to the hadronic ones through the relations 



mx cosh y = zpx cosh yf and mx sinh y = qx sinh 



(6) 
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The inclusive cross section for hadron production can then be written as 



da 



AB^h+X 



da 



AB-*h+X 



r da AB ^ f+x r l 

K(V7s)Eld PT dy f dzD^ h {z,&) 



dq T dy 



dm T dy 



,2 n/ r2 



dp T dy j 



xS(m T - M T (q T , y f , z))5(y - Y(q T , y f , z)) 



(7) 



where the transverse mass and rapidity of the hadron are expressed in terms of the 
integration variables as 



m t(Pt, Vf, z) = (zp T ) 2 + m 2 tanh 2 y f , 



Qt 

Y(pT,yf,z) — aTsmh( — sinhy/). (8) 

TflT 



For the fragmentation functions Df^h we shall use the set KKP [|J which is the latest 
one. The fragmentation scale is //p ~ q?. 

One of the tasks below will be to determine the .fT-factor based on charged-hadron 
production measured in pp(p) collisions at collider energies in the range yfs = 63 . . . 1800 
GeV. The .fT-factor thus effectively accounts for the higher-order effects in the partonic 
cross sections, in the running coupling, in the parton distributions and fragmentation 
functions, and also for the possible intrinsic transverse momentum effects in the parton 
showers both in the initial and in the final state. We shall analyze the yfs dependence 
of K, in particular. 

The parameter p in Eq. ([7]) is the smallest transverse momentum of parton scat- 
terings allowed. For models with semi-hard parton production, this parameter plays a 
key role. We shall also discuss the lower limit for po for each y/s. 

For completeness, we write down the cross section above in an explicit form, which 
can be numerically evaluated in a straightforward manner. Performing the integrations 
over p\ and yj in Eq. (0) gives 



da 



AB^h+X 



dq T dy 



/dz 
— D f ^ h (z,fi 2 F ) 



da AB^f+X 



dp 2 T dy f 



(9) 



f dz f I I 

= K(^ S ).J( mT ,y)J-Jdy 2 ^ — — 

Xlfi/A{Xl, Q 2 ) X 2 fj/ B {x2, Q 2 



(ij)(kl) 



+x 1 f j/A (x 1 , Q 2 ) x 2 fi/B(x 2 , Q 2 ) 



da ij ^ kl da ij ^ k l 
—x {t,u)D k _ h (z,f/ F ) + —j (w,t)A^0,/4) 
dt dt 

da ij ^ kl ~ da ij ^ kl 

-7T (u,i)D k ^ h (z,n 2 F ) + —t (t,u)D l ^ h (z,fjj F ) 

dt dt 



■(10) 



where the partonic cross section is evaluated at 



Pt = —J(rriT,y) } y/ = arsinh(— ^ sinhy) 
z •' q T 



11) 
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and where 



J{m T ,y) = (1 



m 



cosh 2 y ' 



-1/2 



'12) 



The integration region for y 2 , — \og{-\fs/pT — e Vf ) < yi < log(^/s/pT — e Vf ), is over 
the whole phase space, whereas that for z, 



2rriT 



coshy < z < min 



1, —J{m T ,y) 
Po 



(13) 



is affected by the requirement pt > Po- In this way, we shall also be able to consider 
the region q T < p . 

The experiments most often give the cross sections averaged over a pseudorapidity 
interval At/, defined as 



da h AB 



dq^drj 



,, e Ar; At/ J Art dq^drj At/ JArj 



dr]J(m T ,y) 



dn h 
-1 aa AB 



dq\dy ' 



(14) 



where y = arsinh(^ sinliT/) and J(mr,y) 1 = dy/drj, where J(my,y) is the same 
factor as in Eq. flT2l). 



3 The analysis and results for pp(p) 

The approach of collinear factorization and independent fragmentation can be expected 
to work best at the highest values of and at the highest cms-energies. We determine 
K for each based on the data in the large-g^ region by minimizing the x 2 : 

X 2 (N) = E [ \J l J (15) 

where cr^ xp is the measured cross section at a transverse momentum q?i, with a statis- 
tical error A^ xp . The cross sections computed at qxi without the K-factor, are denoted 
by af 1 . The iV data points included in the fit are counted starting from the highest qT 
measured for each y/s. As will be seen, the i^-factor usually depends on N. 

Since in the region of large qx one is sensitive to the (statistical) fluctuations in the 
experimental data, it would be preferable to include as many data points in the fit as 
possible (most preferably all points of course). The previous comparisons [§, |J of the 
pQCD cross sections with the data indicate, however, that the pQCD results tend to 
fall somewhat below the data in the region of qx ~ 3 ... 5 GeV. Also, the pQCD cross 
sections clearly become too large at qr < 1 GeV, which region is sensitive to the cut-off 
parameter pq. The statistical error bars of the data in the few-GeV region are clearly 
smaller than those at large-^T- Thus, if the few-GeV region is included in the fit for 
the determination of the i^-factor, a bad fit for the large-gr region results. This can 
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be seen in Fig. |1] (bottom right), where we plot the i^-factors obtained for the UAl 
MIMI data at = 630 GeV as a function of the smallest transverse momentum q?N 
included in the fit. In the same figure, we also show x 2 /N as a quantitative measure 
of the goodness of the fit. 

To circumvent the problems mentioned above, and to make sure that the i^-factor 
gets determined from the region where the collinearly factorized pQCD approach is 
supposed to work best, we determine the maximum number of data points that can 
be included by requiring a good fit, x 2 (N)/N < 1, at q T > qxN, and fix the ^-factor 
based on these N data points. We estimate the statistical significance of the ^-factor 
at each by computing an error |2T| 



A K {N) 



N ( 3K 

1 1 \da, 



cxp 



(16) 



This gives the error bars shown in Fig. [I], lower right. For the UAl MIMI data we find 
that x 2 (N)/N < 1.0 at q T > 7.988 GeV, with N = 15, and K = 2.19 ± 0.11. 

The qr spectrum of charged hadrons in pp at yfs = 630 GeV, computed from Eq. (TTJ) 
by applying the i^-factor determined above, is shown in Fig. [I] (left panel) together with 
the UAl MIMI data for p + p — > h + X (h = h + + h~). The sum of charged hadrons 
always includes pions, kaons and protons. Both the data and the computation are 
averaged over —3.0 < r\ < 3.0. Contrary to what perhaps could have been expected, the 
agreement with the data are relatively good, taking into account that we discuss particle 
production which extends over 10 orders of magnitude. A more detailed comparison is 
done in the top right panel of Fig. [I], where we plot the data divided by the computed 
cross section, together with (relative) statistical error bars of the data. From this figure 
we can see why large x 2 /N follows if the few- GeV region is included in the fit. 

Fig. [I] (left and top right) also illustrates how we can estimate the lower limit for 
the parameter p : we require that the computation must not overshoot the measured 
cross section in the region q T > p . For this particular data set, we notice that e.g. 
for po = 3 GeV the computed spectrum remains below the measured one until the 
very last data point qx = 0.079 GeV. The value po = 1 GeV in turn would cause an 
overestimate of the cross section already at qr ~ 2 GeV. This search procedure leads 
to a lower limit p = 1.8 . . . 2.0 GeV for = 630 GeV. 

Sensitivity of the computed cross sections to the choice of the fragmentation scale 
Hf and to the different sets of parton distributions is studied in Fig. g for the UAl 
MIMI data. The results are shown for the LO sets CTEQ5 |2|, MRST(c-g) |23|, 
GRV98 Pl| with the scales fiF = Qt and \if = <?t/2. Notice that in this figure we have 
already included the i^-factor which is determined as explained above. The actual 
value of the iT-factor is always correlated with these choices, and for e.g. CTEQ5 and 
/ip = qr/Z, we get K = 1.76 ± 0.05. We notice that qualitatively the behaviour of 
the computed results relative to the data is the same as for fi F = It- Quantitatively, 
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Table 1: The data sets used in the analysis. Systematic errors are quoted when given. The 
if -factors are obtained by requiring a good fit at the high-gr part of the spectra (see the 
text), and the error AK is computed from Eq. (fl^). Systematic error is not included in 
the estimates for AK. Notice that these i^-factors are obtained with the CTEQ5 [^2| set of 
parton distributions, and KKP [|J] set of fragmentation functions and with the scale choices 
fJ-F = Qt, Q = Pt- The values of qxN and N are also shown (see the text). 

the overall agreement between the computed and measured cross sections could be 
somewhat improved by optimizing the choice of the fragmentation scale and to some 
extent also by a choice of the parton distributions (modern sets are to be used, of 
course). For \ip = qx the difference between the best and worst fitting curve is ~ 20% 
and for fi F = gr/2 ~ 30 — 40%. Although not studied here in more detail, there is also 
some freedom for choosing the factorization and renormalization scales, which could 
improve the quality of the overall fit. We do not, however, pursue this study into the 
direction of optimizing the scale choices based on fits to the data but choose to use 
Q = fi = pt, = Qt and the CTEQ5 parton distributions with = 192 MeV in 

what follows. 

The data sets we study are listed in Table 1. The systematic errors given by the 
experiments are quoted. Table 1 also summarizes our results for the fT-factors along 
with their error estimates (systematic errors are not included), and the corresponding 
qxN and N, as well as the obtained values of the lower limits of the cut-off scale p . 

Fig. | presents the pQCD results with the UA1 MIMI data at = 630 GeV 
averaged over various pseudorapidity intervals. We obtained po = 1.9 GeV and K = 
2.19 ± 0.11 using the data averaged over a range \r]\ < 3.0. We use these values to 
decompose the qr spectra into pseudorapidity intervals \r)\ < 0.6, 0.6 < \r]\ < 1.2, 
1.2 < \t]\ < 1.8, 1.8 < \r]\ < 2.4, 2.4 < \rj\ < 3.0. As seen in the figure, these results 
also agree nicely with the data except for 1.8 < |r/| < 2.4. 

Figs. [|and |5] show the computed inclusive charged-hadron cross sections for p + p —>■ 
h + X (h = h+ + h~) at y/s =1800 and 630 GeV together with the CDF data j|. 
Averaging is done over a rapidity interval \y\ < 1.0. For y/s = 1800 GeV (Fig. |], 
bottom right) the distribution of the obtained K(qTN) is nearly flat which indicates a 
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good agreement with the data. Fig. |], top right, shows that this is indeed the case. The 
computed absolute spectrum is shown in the left panel of Fig. |j. The corresponding 
analysis for y/s = 630 GeV is shown in Fig. |^. The agreement with the data seems 
to be again fairly good although the .fT-factor is now fixed only by the last two data 
points. Note that the data only cover a very limited gr-range, 0.425 < qr < 3.5, so in 
practice half of the computed spectrum is controlled by K and the other half by po. 

In Figs. |6|, [7| and |8] we present the computed g^-spectra which correspond to the UAl 
data at y/s = 900, 500 and 200 GeVforp + p -+h+X, where h = (h + + h~)/2. The 
data and the computed spectra are averaged over \rj\ < 2.5. The systematic error for the 
data is ±15% (not shown). The overall agreement with the data for y/s = 900 GeV is 
poor for qr < 5 GeV. A reasonable estimate of po in this case is to take the value which 
at least reproduces the data at smallest values of qr- The result for p = 1.6 GeV with 
the data is shown in Fig. |6| (left panel). The results for y/s = 500 GeV are presented 
in Fig. [7[ The agreement with the data is now reasonable. Also for y/s = 200 GeV, 
shown in Fig. |8|, the agreement with the data is again tolerable. 

The last (and the smallest y/s) dataset studied is the AFS data for p+p — > h + X 
(h = h + + h~) at y/s = 63 GeV, averaged over \y\ < 0.6. The systematic error for 
the data is ±20% (not shown). The results are illustrated in Fig. |[ The data are for 
qr > 2.25 GeV so they cannot be used to determine the value of p . 

The if-factors obtained for different energies are gathered together in Table 1 and 
in Fig. [TIJ The inner error bars correspond to the statistical errors AK computed 
from Eq. fllBD , and the outer error bars stand for the systematic errors alone (when 
available). Except for the UAl dataset at y/s = 900 GeV, the obtained i^-factors 
decrease systematically with increasing y/s. We also observe that a phenomenological 
fit of the form logK = A ± Blog(y/s) can be justified. As the systematic error bars 
are not available for all cases, we do not perform such a fit here but note that based 
on the numbers explicitly given in Table 1, such a fit can be easily performed. 



4 AA-collisions; nuclear effects 

Next we study inclusive charged-hadron production in nuclear collisions. Our emphasis 
is at the highest cms-energies and at the highest transverse momenta. We do not make 
an attempt to invoke any model for the Cronin effect which is observed in pA and AA 
collisions at lower energies P6 |. Different phenomenological approaches to include the 



Cronin effect to the pQCD-motivated framework can be found e.g. in |1C| |Tl| . We do 
not include any modifications of the fragmentation functions, either, which are expected 
to arise in a very dense medium |27], |28|, |29| . Instead, aiming at pQCD cross sections, 
where additional model details are kept in a minimum, we remain within the collinearly 
factorized leading-twist framework where the nuclear effects in hadron production arise 
only from the nuclear modifications of the parton distributions functions. 
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We define the average number density distribution of a parton flavour i in a nucleus 

A as 

Af i/A (x, Q 2 ) = Zff A (x } Q 2 ) + (A - Z)f? A {x, Q 2 ), (17) 

where the nuclear parton distributions (nPDF) are defined in terms of the correspond- 
ing distributions f\ in the free proton as 

f! /A {x,Q*)=Rt{x,Q 2 )h{x,Q 2 ). (18) 

The nuclear modifications Rf depend on the nucleus, the parton flavour, x, and through 
the DGLAP evolution || of ff /A also on the scale Q 2 . For Rf{x,Q 2 ) we use the 



EKS98-parametrization [IH] which is based on a global DGLAP analysis of the nPDF 



31 . The distributions in bound neutrons are obtained through isospin symmetry; 



fu(ti) = fdii) anc ^ fdli) = fu{t)- We ^° no ^ consider the impact parameter dependence 
of the nPDF |32|, |33| here. Clearly, there are two kinds of nuclear effects in the inclusive 
nuclear cross sections to be computed: effects due to the nuclear modifications of parton 
distributions (let us call this shadowing for brevity) and isospin effects. 

The inclusive cross sections at y = in pA and AA collisions at RHIC and LHC 



energies are computed from Eq. (10) with the nPDF defined above. Dividing the 
obtained cross sections by the ones which do not include shadowing but include isospin 
effects, we obtain the ratios shown in Fig. [11] by the dotted and dotted-dashed curves. 
In order to see the magnitude of the isospin effects, we compare the nuclear cross 
sections computed with shadowing and isospin effects against the ones computed for 
pp. The resulting ratios, R p a((1t) and Raa^t) are shown by the dashed and solid lines. 
We can see that the isospin effects remain quite small in all cases. 

The excess in the ratios R p A{qr) and Raa(qt) in Fig. |ll| is caused by antishadowing 
in the nPDF (see Fig. 3 of Ref. |16[). The depletion at small values of qr is due to 



shadowing at small x, and the depletion at larger values of qr is due to the EMC effect 
in the nPDF. The systematics in the location of the excess in g^, its magnitude (the 
height is decreasing with -^/s) is easy to understand by the following argument. 

Approximating the partonic cross section roughly by a power-law, da^ /dpxdy ~ 
CfPr™, and neglecting all masses, Eq. fllCf ) simplifies into 



da NN^h+X ^ £ 

dqrdy ~ f £^ M q% 



E ^ fdzz^Df^q 2 ,), (19) 



For y/s = 200 GeV at p? ~ 10 GeV, we find n ~ 7 and somewhat less for smaller values 
of pt or larger values of y/s. From Fig. [12] for z 6 D(z, g§>) at qr = 10 GeV, we observe 



that the integrand is strongly peaked around Zf = 0.6 . . . 0.8, depending on the flavour 
of the fragmenting parton. Obviously then J dzz n ~ 1 Df^ h (z) ~ z^" 1 D f-*h(z fi <?t) • Con- 
sider only gluons for simplicity (gluons dominate parton production when iCl but 
this is to some extent compensated by the larger fragmentation functions of the quarks). 
We find that da h /dqxdy ~ da 9 /dpxdy where px ~ qr/z with z ss 0.6. 
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At y = 0, the inclusive parton production dominantly probes the region x ~ 2px/ \fs- 
The ratio of shadowed and non-shadowed cross sections in Fig. |Tl| thus behaves as 



^(gT,v^)~[^ A (-^,-)F, (20) 

where, due to the smearing caused by the integration over 2/2, an effective power l<p<2 
appears. Comparison with Fig. 3 of Ref. |31] for R^(x,Q 2 ), shows that this is indeed 
the case: the cross-over points in Raa(<1t) can be traced back to the cross-over points 
of the ratio in Eq. fl2"0|). The decrease in the magnitude of the excess in Raa with 
growing yfs is due to the scale evolution of Rf: at higher scales the antishadowing of 
R^ decreases. Notice also the difference in ^/s for pA and A A at the LHC. 



5 Comparison with the PHENIX data 

As our final task - armed with the yfs dependence of the i^-factors - we compare 
our results with the transverse momentum distributions of charged hadrons in Au+Au 
collisions at ^Js = 130 GeV at RHIC, which have been recently measured by PHENIX 



|r7fl . We consider an average quantity 



dN AA->h+X> , da NN->h+X x x d(J NN^h+X 



(PqTdri I c \ AA ^ ^ d 2 qTdr] ) c ^ binary ^ c a^ N d 2 qrdri ^ ^ 

where Taa is the standard nuclear overlap function, NN refers to nucleon-nucleon 
collisions, and ) c is the average number of binary collisions within the cen- 



trality selection c. In Ref. |]17 |, two centrality classes are given, c\ = . . . 10% and 
c 2 = 60 . . . 80%, along with the estimates (iV binary ) Cl = 905 ±96 and (N hin£LTy ) C2 = 20 ±6 
(resulting from a Glauber analysis) and for the inelastic nucleon-nucleon cross section 
a ia = 40 ±3 mb. We use these values here. To roughly estimate the impact-parameter 
dependence of shadowing, we find an effective nucleus A e g for which the number of bi- 
nary collisions in a central A e ffA e s collision equals (-/Vbmary)c 2 - With Taa(0) ~ A 2 /nR A , 
this leads to A^ ps 163 and A c * s pa 9 for the two centrality classes studied. Shadowing 
is included according to Eqs. ([17]) and ( |18D for A = A eS . Thus the nuclear effects for 
the sample C\{c2) remain within 15(5)% at qx < 10 GeV. 



From the systematics of the -y/i- dependence of the i^-factors in Fig. |10] it is observed 
that a simple power-law interpolation between y/s = 63 GeV and 200 GeV gives a 
fair first estimate of the i^-factor at 130 GeV. This results in K(yfs = 130 GeV) = 
4.27 ±0.18. For the error estimate, we have added the statistical and systematic errors 
of K(tJs = 200 GeV) in quadrature. Applying this i^-factor, we plot the hadron 
spectra shown by the solid curves in Fig. [T3|. For obtaining an estimate of the total 
relative error (shown in the figure with the dashed lines), we add the relative errors of 
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K, (Nftfagjyjc, and a([ N in quadrature. This leads to ±22% and ±36% for the central 
and peripheral cases. 

A very good agreement is found between the computed cross sections and the pe- 
ripheral (c 2 ) data at the highest qx measured. The computed spectrum remains also 
surprisingly close to the measured one at qx ~ 2 GeV which was the trouble region 
in the pp(p) case. The estimated error band now covers the most of expected magni- 
tude of the deviations (see Fig. ||). Within the uncertainties estimated, the peripheral 
Au±Au collisions thus seem to behave effectively as AW collisions. In central collisions, 
however, especially towards the largest transverse momenta, the computed spectrum 
systematically lies above the data. We emphasize that it is in this region that our 
approach should work the best. It should also be noted that if the Cronin effect for 
high-py parton production becomes more important at central Au±Au collisions than 
in the peripheral case, it should increase the computed cross sections. This would in 
turn make the difference between the measured and computed spectra even larger from 



what is seen in Fig. 13. The high-g^ hadron production in central Au±Au collisions 



thus looks quite different from a straightforward extrapolation of the pp case. 



6 Discussion 

In this paper we have studied inclusive charged-hadron production in pp(p) collisions 
at ^/s > 63 GeV in the framework of leading- twist lowest-order perturbative QCD. 
Modern lowest-order sets of parton distributions and fragmentation functions are used. 
Comparison with the AFS |||, UA1 0, UA1 MIMI [|, and CDF [| data is made, and 
the difference between the computed and measured spectra is quantified in terms of 
a factor K = dcr cxp / daf^ Q . We have extracted the K- factors and their a/s- dependence 
from the highest transverse momenta of the measured spectra. A systematic decrease 
of K with growing yfs is found. Also error estimates for the obtained K are given. We 
emphasize that the values of K depend on the PDFs, fragmentation functions and the 
scale choices used. The results of our analysis are useful for further studies of hadron 
spectra, at cms-energies where no pp data are so far available. 

We have also studied the lower limits for partonic transverse momentum exchange, 
the cut-off scale po. A systematic increase of po with growing yfs is found and po = 
1.6 . . .2.2 GeV for i/i = 63 . . . 1800 GeV. It will be interesting to study the relation 
of the scales po obtained here to the saturation scales in AA collisions For the 



models with semihard (perturbative) and soft (nonperturbative) components for par- 
ticle production, our analysis thus suggests a cut-off scale po which increases with \/s. 
The same observation was also made e.g. in a recent analysis for the HIJING model 
Notice that the values of po we obtain are somewhat below those obtained in |T3 



since here we have described all particle production with the perturbative component 
only. 
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Once the i^-factor in each case has been determined, we find the overall agreement 
between the computed and measured spectra in pp at ^/s > 200 GeV relatively good 
(except for the UAl dataset at 900 GeV). In general, the fragmentation functions of 
gluons are not as well constrained by the e + e~ data as those of quarks and antiquarks. 
In the future it will be very interesting to see whether the agreement of the collinearly 
factorized cross sections and the data will be improving with new sets of fragmentation 
functions. It will also be interesting to see how close the results obtained here are with 
the scale-optimized NLO results at various energies. 

The agreement with the pp data at high transverse momenta gives us confidence in 
that a reference cross section can be obtained by making the extrapolation from pp to 
AA collisions. The computed spectra are compared with the recent measurements by 
PHENIX |p.7| . We have shown that in the region = 1...10 GeV for central Au+Au 
collisions at RHIC the antishadowing effects in parton distributions enhance the spectra 
but only by less than 15%. We have also shown that the isospin effects are small. Using 
the estimated K = 4.27 for = 130 GeV, together with the results from a Glauber 
analysis of PHENIX, we have computed the inclusive hadron spectra corresponding 
to the peripheral and central data samples. The results found are very similar to 
those obtained by PHENIX: the peripheral collisions seem to behave practically as 
pp collisions, and the computed spectrum agrees nicely with the data at the largest 
values of qr- In comparison with the data with a 0..10% centrality cut, the computed 
spectrum is found to systematically lie above the measured spectrum at the highest 
transverse momenta, where the emphasis of our approach is. An additional Cronin 
effect, if important at all at this high ^ and large qx, can be expected to enhance 
the computed spectrum and thus enhance the deviation from the data. In order to see 
the effect even more clearly, we are looking forward to more data points at even higher 
values of qx- 



The very dense partonic medium produced (see the initial conditions e.g. in [34 



can be expected to be responsible for the suppression of the high-g^ hadrons p7 |. 
There is an increasing activity to study the fragmentation functions modified by the 
presence of the partonic medium |3R |55 |. Incorporating such modifications into 
the present analysis will be done next [j36) . This procedure focuses on the large-gr 
region, where perturbative methods apply. The other extreme case, emphasizing the 
smaller-g^ region, is hadron production from a fully thermalized system described in 
terms of relativistic hydrodynamics f37|. It will be very interesting to study how these 
two regions merge to form the measured spectrum. 
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Figure 1: Left: Inclusive cross section for charged-hadron production (h = h + + h~) in pp 
collisions at yfs = 630 GeV averaged over \rj\ < 3.0. The LO pQCD prediction with K = 2.19 
and scales Q = pr, = Qt an d po = 1-0,1.9,3.0 GeV is shown by the curves. The data 
shown with the statistical error bars are from UA1 MIMI Top right: The ratios data- 
to-data and data-to-theory as a function of transverse momentum for various pq. Bottom 
right: The mimimized \ 2 (N)/N (solid curve) and the resulting X-factor (with error bars) 
as a function of the smallest transverse momentum, qTN, included in the fit. The X-factor is 
read off from the point where x 2 (N)/N = 1. The errors of the K- factor are computed from 
Eq. (ID. 
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q T (GeV) q T (GeV) 

Figure 2: Left: An example of the dependence of the inclusive pQCD cross sections (as 
in Fig. H|) on the choice of the fragmentation scale fi F . Based on the fit to the large-t/T 
region of the UA1 MIMI data in Fig. we obtain K = 2.19 ± 0.11 for fj, F = q T and 
K = 1.76 ± 0.05 for fi F = q T /2. The scale p has been fixed to 1.9 GeV and CTEQ5 PDF 
are used. Right: As in the top right panel of Fig. |l| but showing the dependence of the 
pQCD results on the PDF set. Three different sets are used, and the scales are fixed to 
Q = p T and p = 1.9 GeV in both panels. For the sets CTQE5L, GRV98 and MRST(c- 
g), we obtain K = 2.19 ± 0.11,2.00 ± 0.10,2.15 ± 0.10 for /j, F = q T (upper panel) and 
K = 1.76 ± 0.05, 1.53 ± 0.07, 1.65 ± 0.08 for /jl f = q T /2 (lower panel). 
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Figure 3: The LO pQCD prediction for the charged-particle spectra in pp collisions at 
y/s = 630 GeV with K = 2.19 and po = 1.9 GeV. The upper left panel corresponds to Fig. |], 
in the other panels the pseudorapidity interval \r]\ < 3.0 has been divided into different 
subintervals. The data are from UAl MIMI, Ref. |§. 
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Figure 4: As Fig. g but for /i = + ft - )/2, y^s = 1800 GeV and |y| < 1.0. The data are 
from CDF gfl. 
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Figure 5: As Fig. | but for yfe = 630 GeV. 
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Figure 6: As Fig. g but for /i = + h _ )/2, = 900 GeV and \rj\ < 2.5. The data are 
from UA1 0. 
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Figure 7: As Fig. | but for ^/i = 500 GeV. 



22 




q T (GeV) q T (GeV) 



Figure 8: As Fig. | but for ^/i = 200 GeV. 
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Figure 9: As Fig. [T] but for h = h + + h , y/s = 63 GeV and \y\ < 0.6. The data are from 
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Figure 10: The obtained if-factors for y/a = 63, 200, 500, 630, 900, 1800 GeV. The inner error 
bars are AK calculated from Eq. ([!(]) and the outer error bars (when shown) are the reported 
systematic errors for each data set (see Table 1). 
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Figure 11: The ratios of transverse momentum spectra of charged hadrons in AA(pA) col- 
lisions relative to pp collisions and AA(pA) collisions without shadowing. At RHIC both 
Au+Au and p+Au are at yfe = 130 GeV and 200 GeV. At the LHC Pb+Pb is at yfa = 5.5 TeV 
and p+Vb at *fs = 8.8 TeV. 
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Figure 12: The integrand of Eq. ( |l9| ) with n = 7. The fragmentation functions are from 
KKP 0. The integrand is peaked at zj. 
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Figure 13: Charged-hadron spectra for centrality classes 0...10% and 60. ..80% in Au+Au 
collisions at ^fs = 130 AGeV at RHIC. The data are from PHENIX @. The solid curves 
are the pQCD results computed with K = 4.27 and A e s = 163 and 9 for the two centrality 
selections. The error bands shown by the dashed lines correspond to an estimated total 
uncertainty of ±22% and ±36%, respectively (see the text). 
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